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Abstract 

The concepts of Wigner time delay and Wigner-Smith matrix allow to characterize temporal aspects of 
a quantum scattering process. The article reviews the statistical properties of the Wigner time delay for 
disordered systems ; the case of disorder in ID with a chiral symmetry is discussed and the relation with 
exponential functionals of the Brownian motion underlined. Another approach for the analysis of time delay 
statistics is the random matrix approach, from which we review few results. As a pratical illustration, we 
briefly outline a theory of non-linear transport and AC transport developed by Biittiker and coworkers, where 
the concept of Wigner-Smith time delay matrix is a central piece allowing to describe screening properties 
in out-of-equilibrium coherent conductors. 
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1. Introduction 

The purpose of this article is to review several 
results linked to a concept which has been very in¬ 
fluential in the work of Markus Biittiker, namely 
the concept of Wigner time delay, and several of its 
extensions - traversal time, Wigner-Smith matrix, 
injectance, etc. In my opinion, this topic allows one 
to have a good flavor of Markus Biittiker’s style as a 
physicist : a combination between formal develop¬ 
ments with mathematical elegance and motivations 
from pratical questions of fundamental condensed 
matter physics. Besides, the choice of this theme 
is also related to a more personal anecdote as it 
was the subject of my first scientific exchange with 
Markus Biittiker, which I rediscovered during the 
preparation of this article, finding in a box a copy 
of an old email, dated the 8th July 1997, addressed 
to me and my PhD advisor Alain Comtet. Here are 
few sentences extracted from this message : “/ no¬ 
ticed yesterday your paper on the cond-mat network 


* Contribution to a special issue Frontiers in quantum 
electronic transport - in memory of Markus Biittiker” 
Email address: christophe.texier@u-psud.fr 
(Christophe Texier) 


m I like to react to two things: (...) In your 
introduction you jump from the work of Fyodorov 
and Sommers immediatly to the work of Brouwer et 
al. Both of these are of course fine works. But the 
work of Brouwer et al. as they make clear in very 
strong terms takes its starting point from my work 
by Gopar and Mello JMi, and that alone should have 
been a good enough reason to not simply leave it 
out. (...) Therefore, it is my hope that you can give 
this work the place it deserves. Sincerely, Markus 
ButtikeF. Of course, he was absolutely right and 
we soon after amended our paper according to his 
remarks. Despite this somewhat awkward begin¬ 
ning, Markus Biittiker hired me as a postdoctoral 
assistant at the University of Geneva two years later 
and I enjoyed very much the pleasant and stimulat¬ 
ing atmosphere of the Physics Department. I espe¬ 
cially appreciated the freedom which Markus gave 
to his postdocs and learnt a lot from our regular 
exchanges. 

The notion of time delay was introduced in the 
fifties by Eisenbud and Wigner in the context of 
the quantum theory of scattering. It allows to 
capture temporal aspects of the scattering pro¬ 
cess [BSl 119311187L HU I122j . The most fundamen- 
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tal aspects of time delays have been reviewed by 
Carvalho and Nussenzveig [63]. Some other articles 
have reviewed more specific aspects : Beenakker 
has considered the case of wave guides in the lo¬ 
calised regime [T3]. In Ref. |107j . Kottos has re¬ 
viewed the random matrix approach for ergodic sys¬ 
tems (chaotic cavities) and the case of (non ergodic) 
disordered systems -diffusive or strongly localised. 
We have emphasized in |54] the connection with 
the theory of exponential functional of the Brown¬ 
ian motion. The purpose of the present article is 
to review few results on time delays from different 
contexts and to emphasize the diversity of the phys¬ 
ical situations where this concept has found some 
applications. 

The outline is the following : in Section we 
discuss several definitions of times characterizing 
the scattering process. Section [^ reviews the sta¬ 
tistical analysis of time delays for disordered sys¬ 
tems. In Section [^ we discuss few aspects of the 
random matrix approach. Section [^introduces sev¬ 
eral generalized concepts (partial DoS, injectance, 
etc.) which will find pratical applications for the 
analysis of non-linear transport (Section]^ and AC 
transport (Section]^. 


2. Wigner time delay and other characteris¬ 
tic times 

2.1. Scattering on the half line 

The most simple situation allowing to introduce 
the concept of time delay is the scattering problem 
on the half line (this also corresponds to project a 
rotational invariant problem in an orbital momen¬ 
tum channel, in higher dimensions). We consider 
the Schrodinger equation 

- 1 p”{x) + V{x)lpeix) = etpe{x) (1) 

for X G K"*" (with f?jlflrn) = 1), describing the 
scattering of a particle by a potential defined on a 
finite interval [0, L] (Fig.[^. At a; = 0 we choose to 
impose a Dirichlet boundary condition, V'e(O) = 0- 
In the “asymptotic” region {x > L), the stationary 



scattering state of energy e = is the superposi¬ 
tion of an incident wave and a reflected 

wave r ■ 


Ipsix) = 


y/2'IThv 




-ik(x-L) _|_ 


( 2 ) 


The normalisation constant, involving the group ve¬ 
locity V = {\/h)de/(ik = 2k, corresponds to asso¬ 
ciate a measure de to the eigenstate. 0 The re¬ 
flection probability amplitude has unit modulus as 
a consequence of current conservation : r = 

The Wigner time delay is defined as the derivative 
of the reflection phase 


Tw(e) = —ihr* 


dr 

de 


^d5r{e) 

de 


(3) 


It measures the time spent by a wave packet of 
energy e in the domain [0, L]. This can be eas¬ 
ily understood by considering the time evolution 
of a wave packet <i)(a:;t) = de $(e) ■!/'e(a;) 
where d>(e) is a narrow function centered around 
an energy eg (from now on, we will set h = 
1 for simplicity). We can split the wave 
packet into an incident part and a reflected 
part : d>(a;;t) = $i„c(a;;t) -b ^ref{x;t). Ne¬ 

glecting the effect of dispersion, ^inc{x;t) = 
(l/V4)()/“dee-i/‘^$(e)e-‘v^("^-^)-'^‘ can be 
rewritten under the form $inc(a:;t) oc f{—x + L — 
vot), where vq is the group velocity at energy eg 
and / the function f(X) « / dK ^{eo+vgK) . 
In the reflected wave packet, the expansion of the 
phase shift produces a shift in time and we obtain 
$ref(a;;t) oc f{x - L - vgt + uorw(eo)) (cf. also 
chapter 10 of |176j b The interpretation of tw as 
the delay in time relies on following the motion of 
a wave packet with sufficiently narrow dispersion 
in energy. As trivial illustrations we can consider 
the case of an impenetrable region when 5^ = tt, 
leading to tw = 0, and a free region {V{x) = 0), 
leading to (5^ = tt -I- 2kL and tw = 2L/r;. 


2.2. The scattering problem in one dimension and 
characteristic times 

The scattering problem on the infinite line il¬ 
lustrates that, although the situation is still sim¬ 
ple, many other characteristic times can be already 
introduced. We consider a plane wave sent from 


^ This choice of normalisation ensures that orthonor¬ 
malisation reads (r/ij Ir/ij,/} = S{e — s') and closure relation 

/delV'eXbsl = 1- 


Figure 1: A plane wave is sent on a potential living on the 
interval [0, L]. 
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as we will explain in the § 2.3 r-w does not have 
a simple interpretation as a wave packet’s delay in 
time. 


Figure 2: The stationary scattering state describes 

the scattering of an incident plane wave from the left, which 
is encoded by a reflection amplitude r and a transmission 
amplitude t. 


—oo (or +oo) on a potential defined on the inter- 
val [-L/2,L/2] (Fig. §. The scattering proper¬ 
ties can be encoded in two stationary scattering 
states (Fig. and 4’e,R{x), controlled by 

two pairs of reflection/transmission amplitudes r, t 
and r', t' characterizing transmission from the left 
and from the right, respectively. A general scatter¬ 
ing state is a linear combination 


^'^(a;) = -b , (4) 


where A = (A^,, A/j) are the amplitudes of the two 
incoming plane waves. They are related to the am¬ 
plitudes B = {Bl, Br) of the outgoing plane waves 
by the 2x2 scattering matrix B = S A : 


S = 



( 5 ) 


As a consequence of current conservation, the scat¬ 
tering matrix is a unitary matrix S £ U{2), 
parametrised by four independent real parameters, 
what is made clear by expressing S in the polar 
representation : 


c _ iVTe-A \ 


( 6 ) 


T = = \t'\^ e [0,1] denotes the transmission 

probability. The three independent phases are a 
global phase $/ = det5, called the ‘^Friedel phase", 
the phase a controlling the left/right asymmetry 
and a magnetic phase x (which can be removed by 
a gauge transformation in the ID case). As the scat¬ 
tering process is characterised by a matrix, there is 
much more freedom to define characteristic times. 


2.2.1. Friedel phase and Wigner time delay 

The most simple generalization of the definition 
@ is 


rw(£) 




1 d^f 

N ~d^ 


( 7 ) 


where N is the number of scattering channels (here 
N = 2). Although this is an important quantity. 


2.2.2. Transmission and reflection group delays 
A plane wave sent on the scattering region from 
the left is splitted into two parts to which may be as¬ 
sociated two phases : a reflection phase Sr = /2+ 

a and a transmission phase St = $y/2-|-x + '^/2- 
Extending the argument exposed in § [2.1 [ would lead 
to introduce three characteristic times (only two for 
a symmetric potential), like the tranmission group 
delay ft = dSt/ds, etc. 


2.2.3. Partial time delays : diagonalisation and 
derivation 

The two scattering states '0 £,q-(x), with a G 
{L, R}, describing a particle incoming from the 
left/right (Fig. [^, may be recombined in order to 
form other basis of solutions. An important case 
corresponds to the two partial scattering states, 
which behave “asymptotically” as[^ 

= [Ca0H{-x)+Da0H{x)]cOs{k\x\+'na) 

ri|>L/2 

( 8 ) 

with o G {1, 2}. 0H(a;) is the Heaviside func¬ 

tion. Decomposing (pg^a{x) over the left/right scat¬ 
tering states, as in Eq. Q, leads to ins] 

det(5 — I 2 ) = 0, i.e. the two partial waves are 

eigenvectors of the scattering matrix with eigenval¬ 
ues We can write : 


where Id is a unitary matrix and {rja} the scattering 
phase shifts. In particular, using the polar decom¬ 
position, one finds the expressions of the two phase 
shifts 

g2ii7i,2 _ gi^//2 _ 7" cos a ± ii/l — (1 — T) cos^ a 

( 10 ) 

leading to 2{r]i +P 2 ) = $/, as it should. Two other 
characteristic times could be introduced, which are 
the two partial time delays 


Ta 


o dr]a 
de ■ 


( 11 ) 


As the partial waves ([^ may be combined in order 
to form wave packets, the partial time delays can 
be interpreted as the delay of such wave packets. 


2 For a symmetric potential V(x) = V{—x), the two par¬ 
tial waves are the symmetric and antisymmetric solutions 
(i.e. Cl = Di and C 2 = —D 2 ). 
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2.2.4- Proper time delays : derivation and diago- 
nalisation 

The determination of the partial time delays in¬ 
volves a diagonalisation of the scattering matrix 
and a derivation of its eigenvalues. The converse of 
these two operations leads to introduce the proper 
time delays {t^}, which are eigenvalues of the 
Wigner-Smith time delay matrix 

8 = (o ") 


where V is a unitary matrix, which differs from 
U involved in ([^ in general (see Appendix B). 
Hermiticity Q = follows from unitarity of S. 
Proper time delays and partial time delays have an 
important difference : being derivatives of the S- 
matrix eigenvalues, the partial time delays are in¬ 
trinsic properties of the scattering process, whereas 
the proper time delays depend on the particular 
choice of basis in which the S matrix is expressed 
(see Appendix B). Although they do not coincide 
in general, they satisfy the sum rule 


d N N 

- i— Indet 5 = fg = Tg = A^Tw 

a—1 a—1 


(13) 


where N is the number of scattering channels (here 
N = 2). 


2.2.5. Other characteristic times 

The time delay interpretation of the various char¬ 
acteristic times deduced from scattering phases may 
lead to paradoxal conclusions in the presence of tun¬ 
neling barriers, like superluminal propagation. For 
this reason “clock approaches” have been proposed 
in order to provide more satisfactory answers to the 
question “how much time needs a wave packet to 
travel in a given region ?”. We refer to the review 
articles [Ml 11161 1551 1115L |32J |S3] and Biittiker’s 
contributions 1551H5] . 


2.3. Krein-Friedel relation and Virial expansion 

An important aspect behind the notion of time 
delay is its link with the spectral properties of open 
systems with continuous spectra, which can be un¬ 
derstood as follows. Let us come back to the simple 
situation of scattering on the semi-infinite line for 


^ The introduction of Ref. m, which introduces a so¬ 
journ time from wave attenuation, provides a good overview. 


simplicity (Fig.[^. The scattering states ([^ satisfy¬ 
ing the Dirichlet boundary condition ^/’e(O) = 0 may 
be used as a basis in order to consider the spectral 
problem on the finite interval [0, L], which involves 
a second boundary condition at x = L. We choose 
a Neumann boundary condition 'fi'fiL) = 0 for in¬ 
stance. The eigenenergies {e„} are solutions of the 
quantisation equation Sr{en) = 2mr for n integer. 
When two successive energies are sufficiently close, 
we can expand the relation Sr{en+i) — dr(Sn) = Stt 
as S^(en)5£n — 27r, where = £n+i — £n is the 
level spacing. We deduce the expression of the den¬ 
sity of states i^(e) ~ 1/fen — Tw(s)/(27rh). A 
more precise connection between spectral and scat¬ 
tering properties was obtained by Friedel m and 
Smith |166j . who established the relation 

^ dx \ fie{x)\^ = ^ (^Tw(e) + (14) 


(see also |179j l. The left hand side of the equation is 
simply related to the local density of states (DoS) 

y{x-,£) = {x\S{£ - H)\x) = . (15) 


Thus, Eq. (14) shows that the Wigner time delay 
provides a measure of the DoS of the interval [0, L], 
and may be written in the more general form 


.-(£) = / d.„(*;s) = TTr|stg + 2^}, 

( 16 ) 

SO that it also applies to more general situations 
such as the one of Fig. or even complex struc¬ 
tures like metric graphs |175Lll78lll80j . Such equal¬ 
ities are known as Krein-Friedel relations _piiina 
iTTl |3l] (or Birman-Krein formula [22] )• |j They 
have a long history and first appeared in the virial 
expansion of the equation of state of real gases 

p = nT[l + B 2 {T) n + Bsir) n^H-], where p is the 

pressure, n the density and T the temperature. As 
shown by Beth and Uhlenbeck |188L [20] , the sec¬ 
ond virial coefficient can be written as B 2 {T) = 


^ The choice of the normalisation for the scattering state 
§ is important. 

^ The Krein-Friedel relation is sometimes written in terms 
of the spectral shift function A(£) = f^^dx[u{x;e) — 
i'o{x;£)^, which measures the variation of the DoS due to 
the introduction of the scattering potential, where i'q{x]£) 
is the free local DoS. The precise relationship between the 
two formulations has been discussed in Ref. [175| for the ID 
situation. When the scattering matrix is defined like in 
we have A(e:) = (2i7r) ^^^[indettS — 2iy^L]. 
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—2 ®/^Ai^(±l+16Zint) for bosons (+) and fermions 
(—), where j (jnT) is the thermal 

length and m the mass of the particles. The first 
contribution to B 2 encodes the quantum correla¬ 
tions arising from the symmetrisation postulate, 
whereas the second term describes the correlations 
due to the interaction, being the partition func¬ 
tion of the two body problem (in the relative coor¬ 
dinates) 

= 5](2£ + 1) / —n{e) , (17) 

e Jo 

where summation runs over orbital momentum. 
T( = 2dr]e/de is the time delay related to the scat¬ 
tering phase shift characterizing the scattering 
in the channel of orbital momentum i (cf. § 77 
of |114j l ; in the presence of bound states, the 
partition function receives additional contributions. 
This approach was later made more systematic by 
Dashen, Ma and Bernstein [S2]. Note that a sim¬ 
ilar analysis of real gases based on time delay can 
be developed within a purely classical frame, as ex¬ 
plained by Ma m- 


2.4- Few remarks 

Characterisation of the spectral properties of 
open systems has important applications as it al¬ 
lows to express several physical observables in 
terms of scattering properties at the heart of the 
Landauer-Biittiker description of quantum trans¬ 
port. For example, (16) allows to measure the 
charge inside a conductor in terms of its scatter¬ 
ing matrix, an idea which turns out to be central in 
many developments of Biittiker involving screening, 
which are reviewed below. Several remarks : 


• The Krein-Friedel relation was studied in met¬ 
ric graphs in |175L 1178) . In such systems with 
non trivial topology, the system might pos¬ 
sess so-called “bound states in the continuum” 
(BIC), i.e. a discrete spectrum superimposed 
onto the continuum spectrum (this occurs 
when symmetries allow some bound states to 
remain uncoupled to the continuum) |192LI171] . 
In this case, the Krein-Friedel relation only 
characterises the continuous part of the spec¬ 
trum |178M180j (see |175[I178) for illustrations). 


® Note that 0 assumes that the scattering phase shifts 
are measured with respect to the free eigenstates, ijjiir) ~ 
as it is done usually in scattering 

theory, instead of 0 


A related problem concerns the role of the 
transmission phase. As it is clear from the po¬ 
lar representation (§, the Friedel phase 
which enters the relation (16) as det S = , 

is related to the transmission phase 6t = {^f + 
7r)/2 in ID, as the magnetic field can always 
be removed by a gauge transformation on the 
infinite line. Thus the DoS can be as well re¬ 
lated to St in ID m- If the system has a 
complex topology, the transmission may van¬ 
ish which causes transmission phase jumps. As 
a consequence of these phase jumps, the Friedel 
phase and the transmission phase differ ; the 
density of states is then only related to the 
Friedel phase, as discussed by Taniguchi and 
Biittiker m- This question has been exam¬ 
ined numerically in quantum dots by Levy Yey- 
ati and Biittiker 11191. 


• In Ref. |168j . the Krein-Friedel relation was 
considered in the quasi-one-dimensional situ¬ 
ation (wave guide) which is relevant in order 
to describe the electric contacts of mesoscopic 
structures. 


3. Disordered systems (localised regime) 


The review has underlined the close relation¬ 
ship between Wigner time delay for ID disordered 
systems and exponential functionals of the Brow¬ 
nian motion, which have been widely studied in 
the mathematical literature |195) (see also [SS] for 
a brief review). Let us first recall few properties 
which will be useful for the following. 


3.1. Exponential functionals of the Brownian mo¬ 
tion 


We introduce the random variable 

, (18) 
Jo 

where B(x) is a normalised Brownian motion start¬ 
ing from B{0) = 0 (Wiener process). The distribu- 
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tion of (181 was found in [1361 155] (see also [SS]) : 




(2Zy+i^ ^ r(l + a-n) 

' ' 0^n<iJ./2 1 ^ ’ 


X e 
= -1/(42) 


+ 


— 2Xn{^—n) r ^—2n f _|_ 1 

^ ’ ” \2Zj 

ds s sinh(7rs) 


r 


IS — /i 
2 


27r2(2Z)(i+M)/2 7o 

X f^(i+M)/2,W2 , (19) 

where L'^{x) is a Laguerre polynomial and W^^^(z) 
the Whittaker function [5^. Eq. (191 shows that 
the random variable admits the limit law 

1 




, e for > 0 . 


2i^r(/r) Z^+^^ 

( 20 ) 

For finite X, the functional is characterised by ex¬ 
ponential moments 

{(ZW) ~ . (21) 

^ ^ ' ' x^oo r(2n —/r) 

The precise expression of the moments can be found 
in Ref. |143j for /i = 0 and Ref. |136j in the general 
case. Q 


3.2. Disorder on the half line 
The interest for time delays in disordered systems 
has started with the work of Sulem and cowork¬ 
ers [79] on stochastic resonances (see also Ezl). 
The idea that localised states could produce sharp 
resonances was later developed by Azbel [T2]. The 
first statistical analysis of the Wigner time delay in 
a one-dimensional situation was provided by Jayan- 
navar, Vijayagovindan and Kumar [104] . These 
authors have studied the time delay distribution 
Pl{x) for the Schrodinger equation 0 on the half 
line (Fig. when V{x) is a Gaussian white noise. 
They have identified the existence of a limit law 
with power law tail Poo{t) ^ for an infinitely 
long disordered region, L —)■ oo, although their 
expression of Poo(t) was partly incorrect (with a 
non vanishing distribution as r —>■ 0). Similar 
conclusions were obtained in Ref. m when the 
potential V(x) is the integral of a white noise, 
i.e. a Brownian motion. As it was recognised 


^ Note however that Eq. 4.4 of Ref. 11361 for /r = 0 con¬ 
tains a misprint : should be replaced by Cf~^. 


later [iss], the incorrect non vanishing Poo(O) of 
Refs. [10411105] was due to some inappropriate av¬ 
eraging over the fast phase variable in the weak dis¬ 
order regime. Short after, Heinrichs identified that 
the exponential moments, (r") ^ exp[2n(n—1)L/^] 
where ^ is the localisation length, are character¬ 
istic of a log-normal tail Pl{t) when the length 
L is finite unoj. Inspired by the work of Faris 
and Tsay [72], we established in Ref. m the con¬ 
nection between the Wigner time delay and expo¬ 
nential functionals of the Brownian motion. We 
analysed two different disordered models exhibit¬ 
ing qualitatively different spectral and localisation 
properties, namely the Schrodinger Hamiltonian 
p[ — —dl + V(x) and the supersymmetric Hamil¬ 
tonian P[ = —df + m{x)‘^ -I- m^(x) [i.e. the square 
of the Dirac Hamiltonian (29)], where the potential 
V(x) and the mass m(x) are Gaussian white noises. 
This analysis suggested the universality of the sta¬ 
tistical properties, which was later demonstrated in 
Ref. m- This can be understood from the integral 
representation ( |14| ) : in the weak disorder regime we 
may write rw — ‘2,tt /q^ da; |'0e(a;)p. Using that the 
logarithm of the envelope of the wave function is 
a Brownian motion with drift, we have related the 
time delay to the functional (181 


(law) (1) 

rw = 2rjZ^/^ 


( 22 ) 


where ^ is the localisation length]^ and = ^/v = 
f/{2.k) is the time needed for a ballistic motion on 
the localisation length The value of the drift 
/r = 1 originates from the equality 

(ln|V’£(a;)|) ~ Var(ln|^/i£(a:)|) for a; > ^ , (23) 


valid in the weak disorder regime [TO]. Eq. ( [2^ 
is known as “single parameter scaling'’ [7| 153] . 
as the full distribution of In !■!/>£ (a;) | (or of the 
conductance) is characterised by a unique length 
scale, the localisation length defined by 1/^ = 
lima;_).oo(l/x) In \'ipg(x)\. This question has been re¬ 
discussed more recently for ID disordered systems 
The universality of the Wigner 


^ An equality in law relates two quantities with the 
same statistical properties. For example, the well-known 
scaling properties of the Brownian motion may be conve¬ 
niently written B(Xx) where B(x) is a Brow¬ 

nian motion. 

® For the model Q with {y(x)y(a:^)} = aS(x — x'), the 
localisation length is ^ ~ Se/a for weak disorder e. 
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time delay statistical properties are thus under¬ 
stood as a direct consequence of the universality 
of the localisation properties in ID in the weak dis¬ 
order regime. E3 In other terms, Eq. ( [2^ means 
that the distribution of the Wigner time delay is, 
in the universal (weak disorder) regime, 


If the length L of the disordered region goes to in¬ 
finity we get the limit law 117411179j 

Poo{T) = ^e-^^/^ . (25) 

The limits L —>■ oo and r —>■ oo do not com¬ 
mute : for finite L the tail of the distribution is log¬ 
normal |100Lll79j . Pl{t) ^ exp [—(^/8L) \'a^(T/T^)\ 
for T —>■ oo. The leading term of the moments was 
also given in Ref. m ■■ 


(t) = 2vL 


for n > 1 


(26) 

(27) 


confirming the main exponential behaviour found 
earlier by Heinrichs m and providing the pre¬ 
exponential factor (a more precise expression of 
the moments was given in Ref. |179j b Note that 
the log-normal distribution is characterised by mo¬ 
ments of the form exp[cn^] ; the fact that the mo¬ 
ments \27\ increase with a lower rate with n, as 
~ exp[2n(n— 1)L/^], reflects the existence of a limit 
law with tail (see § 3.1). The importance of 
Azbel resonances was emphasized as it allowed to 
recover the moments by some heuristic argument 
m- A lattice model was studied in Ref. |145| . 
leading to the same conclusions. 


3.3. Disorder in ID with a chiral symmetry 

The localisation properties of disordered systems 
are mainly controlled by their dimensionality [T] 
and symmetries. A first classification of symme¬ 
tries was given by Wigner and Dyson in the context 
of random matrix theory, depending on the pres¬ 
ence of time reversal symmetry and spin rotational 
symmetry, leading to the orthogonal (/3 = 1), uni¬ 
tary {(3 = 2) and symplectic classes (/3 = 4) |129| . 


The time delay distribution was also studied in the non 
universal (strong disorder) regime in El] for a particular 
model. 


These three symmetry classes were later completed 
by others, identified according to the presence or 
not of two other types of discrete symmetries which 
might occur in condensed matter physics : the 
chiral (or sublattice) symmetry and the particle- 
hole symmetry. These has led to the classifica¬ 
tion in terms of ten symmetry classes [naiaiTi]. 
A one-dimensional disordered model with a chi¬ 
ral symmetry which has attracted a lot of atten¬ 
tion is the Dirac equation with a random mass (see 
PSI1551 1181] for reviews) 


with 




T-Ld = io'29a; -I- (Ti m{x) , 


(28) 


(29) 


where (Ji are the Pauli matrices. A chiral symme¬ 
try is the anticommutation of the Hamiltonian with 
a unitary operator, here = —Pd- The 

Wigner time delay distribution for this model was 
obtained in Ref. m in the “universal” (weak dis¬ 
order) regime, leading to the same conclusions as 
for the Schrddinger equation, Eqs. ( |22p4t , provided 
that the localisation length is modified, ^ ~ 

2/g for £ :§> g (in the Dirac equation, the velocity 
is equal to unity). The energy e can be viewed as a 
parameter which tunes the chiral symmetry break¬ 
ing : at the symmetry point e = 0, where chiral 
symmetry holds, the Wigner time delay presents 
different properties as we now discuss. 

The scattering problem on the half line for the 
Dirac equation is settled as follows : we consider 
again the situation of Fig. where the mass is non 
zero on the interval [0, L\. In the free region we can 
decompose the spinor of energy e > 0 as the com¬ 
bination of an incoming plane wave and a reflected 
plane wave : 


is(x — L)+iSr.{e) 


(30) 

Following [5711174] we can parametrize the spinor on 
[0,1/] as 'k(a:) = p(a:) (sin0(a;), —cos0(a::)), which 
leads to the equation for the phase 


d6{x) 

dx 


= e + m{x) sin [20(a;)] . 


(31) 


The case of higher dimensions is also very important : 
besides that it appears in several contexts of condensed mat¬ 
ter physics (graphene, etc), it is of fundamental interest as 
it was shown that all symmetry classes can be represented 
within the 2D Dirac equation with randomness m- 
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The matching with (30) shows that the phase shift 


is given by Sr = 20 (L). We choose the initial con¬ 
dition 0(0) = 0 or 7r/2, which corresponds to an in¬ 
finite mass m(x) = -l-oo or m(x) = —oo for a: < 0, 
respectively (these boundary conditions confine the 
particle on K_|_ ; they are the only ones which do 
not break the chiral symmetry m)- Following 
m we introduce the variable Z{x) = 2d9{x)/ds 
which obeys 


dZ{x) 

dx 


= 2 -b 2m{x) Z{x) cos [26(x)] (32) 


with initial condition Z(0) = 0 and provides the 
value of the Wigner time delay rw = Z{L). 

The analysis of the symmetry point e = 0 is easy : 
the phase remains locked at 0 or 7r/2 (corresponding 
to m{x) = ±oo for a; < 0), therefore 


tw = 2 da:e±2/i'd^'™(^') 

Jo 


at e = 0 . (33) 


For the sake of concreteness, we now consider 
the situation where the mass is a Gaussian white 
noise {m{x)m{x'))^ = gS{x — x') with mean value 
{in(x)) = ng. The representation (331 makes clear 
the identity 


Although the distribution has no limit law, the 
large L limit presents the power law behaviour m 


^ L->-oo Y^27rgLr 


-i/ls-r) 


(37) 


cut off by the log-normal tail for r —>■ oo. The 
absence of a limit law may be associated with the 
delocalisation of the model for e = 0, as all time 
scales extracted from the distribution increase with 
L (see [MIISHIIIHI] for reviews). 


3.3.£. The case {m{x)) = gg ^ 0 : time delay as a 
probe for zero mode 

We now analyse the case {m{x)) = gg ^ 0 which 
was not considered in the literature. Although the 
spectral properties of the model are invariant un¬ 
der the change of the sign of the mass, this is not 
the case for the scattering properties : neither for 
the phase distribution [S2] nor for the Wigner time 
delay distribution, as we demonstrate here. The 
distribution (t) has an interesting property : the 
existence of a limit law is correlated with the choice 
of the boundary condition at a; = 0 and the sign of 
the average mass {m{x)) = gg. If we choose the 
boundary condition corresponding to m(x) = ±oo 
for X < 0, we obtain the limit law only for =F/r > 0 


(l™) z 

at e = 0 , 

i.e. the distribution isF^ 



(f). 


(34) 


(35) 


where the dimensionless function was defined 
above, Eq. (19). 


3.3.1. The critical case (to(x)) = 0 

The distribution for g = 0 was obtained in m 
(although the distribution has been already deter¬ 
mined in another context m, what was used 
in |174] 1. It is characterised by a log-normal tail 
Pl{t) ^ exp [ — l/(8gL) In^(r)] and exponential 
moments 


It is interesting to compare | |35| | for e = 0 with the high 
energy (e g) distribution ( |24[ l obtained in the absence of 
a chiral symmetry, which takes the explicit form Pl{p) = 

( 9 / 4 ) V’h^/2 (9r/4). 


Poo{t) 


g p-^Kar) 

F(|^|)(gr)i+IH 


(38) 


Conversely, for ±g > 0, there is no limit law as 
L ^ oo. We can correlate the existence of the limit 
law with the presence of a chiral zero mode located 
at the boundary x = 0, what occurs when the sign 
of the mass changes at the boundary : 


m(x) 
for X < 0 

(to(x)) 
for X > 0 

zero 

mode 

limit 

law Poo{t) 

— OO 

> 0 

yes 

yes 

—00 

< 0 

no 

no 

-1-00 

> 0 

no 

no 

-1-00 

< 0 

yes 

yes 


(the existence of zero modes in the multichannel 
Dirac equation with random mass was recently 
studied in Ref. [S2], where their topological nature 
was underlined). 


3 . 4 . Related questions and final remarks 

• An alternative derivation of the stationary dis¬ 
tribution (25) was proposed in Ref. m, by 
using the relation between the time delay and 
the reflection coefficient in the presence of a 


























constant imaginary component in the potential 
(wave amplification), i.e. making use of ana¬ 
lytic properties of the scattering matrix (see 
also [TTlISn] '). 

• A discrete Anderson model with Cauchy dis¬ 
order (Lloyd model |121j ') was studied in 
Ref. [53], for which the relation between the 
two cumulants of ln|'0£(a;)| presents an extra 
factor 2 (i.e. drift /r = 1/2), compared with 
the usual form for the single parameter scal¬ 
ing, Eq. @.0 Some numerical investigations 
have shown that the distribution of the Wigner 
time delay is however still characterised by a 
power law tail (and not r“^/^). It would 
thus be interesting to clarify this observation in 
connection with the discussion of the section. 

• Correlations in energy {T^f/{e)Ty/{e')) were 
analysed by Titov and Fyodorov |185j . 

• Time delay and resonance width are closely re¬ 
lated (see the review |107j L The resonance 
width distribution was studied in [1131196j . 

• We have focused the discussion on the Wigner 
time delay characterizing the reflection prob¬ 
lem on the half line. The time delay describ¬ 
ing the transmitted wave through a disordered 
medium was analysed in [25] . 

3.5. Beyond ID 

3.5.1. Multichannel disordered wires 

The extension to the quasi-lD situation (a wave 
guide with N channels) has been considered by sev¬ 
eral authors (e.g. see |138l 1190] L Using the re¬ 
lation between the time delays and the reflection 
probabilities in an absorbing or amplifying medium 
(i.e. the analyticity of the 5-matrix), Beenakker 
and Brouwer obtained the distribution of the proper 
time delays describing the scattering on a semi¬ 
infinite disordered wave guide [Hdl] : 

'Pi'ji,--- ,1 n) oc ]J|7z 

i<j k 

where jk = Tsl'^k is the inverse of the proper time 
and Tg a characteristic time related to the disorder 
strength. The result relies on the isotropy assump¬ 
tion (among the channels), i.e. describes the weakly 


disordered (quasi-ID diffusive) regime ; in this case 
the localisation length increases with the number of 
channels as [13] ^(q^asi id) _ ^(id) |^1-)-^(]V_ i)/2] . 
The distribution ( [39| is a particular instance of 
the Laguerre ensemble of random matrix theory 
(RMT). For one channel N = 1, one recovers the 
limit law (251 by setting I3tsI2 = t^. 

Starting from (391, the distribution of the Wigner 
time delay r-w = (1/-^) obtained in 

Ref. [Ml : 




Ca 


X exp ■ 


27/3 


64 


1 - 


(40) 


/3\ 9(2-v^) 


4r 


where /3 is the Dyson index for orthogonal (/3 = 1) 
or unitary (/3 = 2) symmetry classes. Cp is a 
normalisation constant. The distribution presents 
the same power law tail as in the strictly one¬ 
dimensional case ^ This can be 

understood from the fact that the physics at large 
time (i.e. large scale) is expected to be dominated 
by a single channel (the less localised one). 

Denoting by vl{£f) the DoS of the multichan¬ 
nel disordered wire of length L, (40) can be inter¬ 


preted as the limit law for tw — 2tt VL{eF)/N in 
the L —7> oo limit. The fact that all moments of 
Tw are infinite in this case is thus simply related to 
the divergence of the DoS when L —)■ oo, as in the 
strictly ID situation. 


3.5.2. Higher dimensions 

• Some measurements of time delays were 
performed with electromagnetic waves in 

Ref. [T5511551155]. 

• The case of higher dimensions was considered 
by Ossipov and Fyodorov |144j who analysed 
a two dimensional situation in the diffusive 
regime. The localised regime was studied by 
numerical simulations in Ref. |194j . 

• A study of time delay at a critical point like 
the metal-insulator transition was performed 
in Ref. [TOHl[80l[l44] (see the review [T7)7] L 


4. Random matrix approach for quantum 
dots 


Localisation for potentials with power law distribution 
was studied in 12111651 Il21l 11861 . 


In systems with ergodic properties, like a chaotic 
cavity with narrow contacts (Figs. and § , or a 
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weakly disordered cavity in the ergodic regime, it 
is natural to make a maximum entropy assump¬ 
tion |132[ [78l [CT I130L I131j leading to postulate that 
the scattering matrix is uniformly distributed over 
the unitary group (for perfect couplings at the con¬ 
tacts). Although this approach has led to many 
successes in the description of several properties 
of coherent conductors (conductance, shot noise,...) 
mm^ it does not provide any information about 
the energy dependence of the scattering matrix, 
which is probed by the Wigner-Smith matrix (12). 
Two approaches were proposed to develop a ran¬ 
dom matrix description for the energy dependence 
of the scattering matrix : 

• the “Hamiltonian Approach” (HA) of chaotic 
scattering pushed forward by Fyodorov, Savin, 
Sommers and coworkers, 

• the “Alternative Stochastic Approach” (ASA), pi¬ 
oneered by the work of Brouwer and Biittiker [28] 
and mostly developed by Brouwer, Beenakker and 
coworkers. 

One of the first result within RMT was the cal¬ 
culation of the two point correlation function of the 
time delay m, Eq. (ro^. The analysis the distri¬ 
butions was considered soon after by Fyodorov and 
Sommers [531 [53] who derived the marginal law for 
partial time delays = {1/N) Y' { S{t — Ta)) 

in the unitary class (/3 = 2), Eq. (|49|rp^ 

The distribution of the Wigner time delay, 
( 1 /^)!]*^)), is of great inter¬ 
est. We also discuss below the marginal law for the 
proper times = (l/N) - (S(t - n)). 


4-1- Wigner time delay for the Wigner-Dyson en¬ 
sembles 

4 . 1 . 1 . Distribution and moments 
Proper time delays {Ti} —. Gopar, Mello and 
Biittiker obtained the Wigner time delay distribu¬ 
tion in the one channel case for the three symmetry 
classes (90] : 

(obviously = w[^\t)). The 

proper times are measured in unit of the Heisenberg 


Note that a similar route was followed in Ref. Cm]. 
however the result of the reference is partly incorrect as ex¬ 
plained in Ell- 


time th = 27rfi./A, where A is the mean level spac¬ 
ing of the cavity (of more correctly the “mean res¬ 
onance spacing” as we deal with a scattering prob¬ 
lem). 

An important advance was the work of 
Beenakker, Brouwer and Frahm who were able to 
obtain the joint distribution of the proper time de¬ 
lays within ASA in Ref. |29l |30| . They established 
the relationship with the Laguerre ensemble of ran¬ 
dom matrices : the inverse of the Wigner-Smith 
matrix is a Wishart matrix Q~^ = where X 

is a N X {2N -1-1- 2//3) matrix with identical and 
independent Gaussian elements (this interpretation 
only holds for /? = 1, 2). The joint distribution of 
the inverse of the proper times, 7 ^ = th/t^ being 
given by [inil3D] 


,7Ar) . 

i<j k 

(42) 

Eor N = 1, the Gamma-law 7 ^( 71 ) oc 
corresponds to Gopar, Mello and Biittiker’s result, 
Eq. (41). 

In the large X limit, the marginal law for proper 
time delays takes the form 


Np{Nt) 

N—>-oo 


where [3nl[T82][g 


P{y) = ^;^\/(a;+ -?/)(?/-a;-) 


(43) 


(44) 


with x± = (v^ ± 1 )^. The marginal law is not 
strictly zero out of the interval [x-/N,x+/N] but 
presents finite N corrections with asymptotics 0 

(45) 

~ . (46) 

r—>-oo 


The distribution is related to the Marcenko-Pastur 
law by p(y) = pyivD/v), where PMp(a:) = 

— x){x — X-) is the distribution of the 
rescaled eigenvalues 'y = N x for N x N matrices in the 
Laguerre ensemble, Eq. ( |42[ |. 

In the general case, the Laguerre ensemble describes 
N X N random Hermitician matrices T with positive eigen¬ 
values distributed according to P(r) oc (det exp [ — 

(/3/2) Tr {r} ] (Eq. \A2\ corresponds to a = 1 and l |39[ l to a = 
0). The Marcenko-Pastur law 'W]s[{'y) ~ (1/-N) pMP(q/N") 
for the density of eigenvalues has support [Nx—,Nx^] 
where x± = (a/T+Tk di l)^. Finite N corrections to the 
Marcenko-Pastur law for 7 G [0, A^a;_] U [N’a 7 +,oo[ were 
obtained by Forrester in Ref. ESj. Asymptotic expansions 
of Forrester’s result are iCjv( 7 ) ~ ryOL^N /2 ^ ^ q 

for 7 —>■ 00 . 
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We can deduce the moments : 

{n) = ^ and Var(T,)ci^ V/3. (47) 

The approach of Ref. m provides the precise ex¬ 
pression of the variance m- 


• In the unitary case, a general formula for the 
marginal law for partial time delays w)^ (r) for 
N arbitrary (fixed) couplings between lead and 
dot was also obtained in Ref. [51], which re¬ 
duces to (49) when all channels are perfectly 
coupled. 


= iV[/?(^-l) + 2]+2 
^ N^{N + 1){PN -2) ■ 


(48) 



y 


Figure 3: Marginal law for the proper time delay. 

Remarks : 

• The role of non ideal couplings was also con¬ 
sidered in Ref. where a generalisation of 

the Marcenko-Pastur law was derived (density 
of eigenvalues in the Laguerre ensemble). 

• Joint moments of proper times were studied in 

m- 


Partial time delays {fa} —. The marginal law for 
partial times was obtained for /3 = 2 in [521 [M] and 
for /3 = 1 in [52]. The distribution V /3 is given in 
Ref. [ISg : 


w 


^ (/ 3 / 2 )^^/^ ^_2-PN/2 ^-f3/(2r) ^ 


Nr{pN/2) 


(49) 


The distribution matches with (41) for N = 1, 
as it should. We deduce the moments (fa) = 
(/?/2)"r(l -n + /3iV/2)/r(l + /3iV/2), i.e. 


(Ta) = 


1 

N 


and V<u'(f.) = (50) 


Thus the fluctuations of the partial time delays are 
reduced compared to those of the proper time de¬ 
lays : Var(fa) ~ 2 Var(Tj)/(/3iV). 

Remarks : 


• The crossower between orthogonal (/3 = 1) and 
unitary (/3 = 2) classes was studied in Ref. [52] 
with the effect of non perfect coupling : in the 
weak coupling limit, when g — 2/T —1^1 
where T is the transmission probability, the 
marginal distribution was shown to present the 
tail ~ intermediate 

range g~^ t g, whereas the far tail is not 
affected ■u;^^(r) ~ -j--' 2 - 0 n /2 


Wigner time delay rw —. Despite the joint dis¬ 
tribution for the proper times was explicitly 
known [23], the statistical properties of their sum, 
the Wigner time delay rw = = 

(1/A^)^jjfa, remained unknown for a while. The 
distribution for N = 2 channels was deduced in 
Ref. [HO] : 


/?^/^+^r(3(/3 + l)/ 2 ) 

2 ^ ’ r(/3-k i)r(3/3-k 2) 

XT-3('5+i)t/(^^,2(^+l);/3/r) 


(51) 


o-Ph 


where U (a, 5; z) is a Kummer function (confluent 
hypergeometric function) [2]. A method providing 
a systematic determination of the cumulants was 
proposed by Mezzadri and Simm in Ref. |135j . The 
authors gave explicitly the first four cumulants 


{tw) 





m 

N 

_ 

iV2(A^-k l)(iV/3 - 2) 


(52) 

(53) 


_96'Tj}_ 

+ 2){N + 1){NP - 2){NP - 4) 

(54) 


96(53Af^-68Af-156)TH 

N‘^{N+3)iN+2){N+iy(N-2y(N-A)(N-6) 

for /3 = 1 

12(53Af^-77)T| 

< N‘^{N+3){N+2)(N+l)^(N-l)'^{N-2)(N-3) 

for /3 = 2 


12(537V^-|-34A)'-39)th 

N‘i{N+3)(N+2){N+iy{2N-iy{N-l){2N-3) 

for /3 = 4 


(55) 
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The analysis of the variance is instructive : we 
can write 


Var(A^ Tw) = Var( ^ n) ~ Var(Ti) 

i 


= Var(^fa) ~ 2A^Var(fa) . 

a 


(56) 

(57) 


Writing Var(A^ r^) = iVVar(ri) + N{N — 

1) Gov{Ti,Tj) and using ( [48| ) provides the covari¬ 
ance Cov(Ti,Tj) = —1/[N^N + 1)], showing that 
the proper times are anti-correlated. We have 


Cov(Ti,Tj) fiN— 2 1 

^Var(Ti) Var(Tj) iV[/3(iV — 1)-|-2] -|-2 N 

(58) 

A similar argument gives the covariance for partial 
times, Cov{fa,h) = 2/[N'^{N l){N/3 - 2)], re¬ 

sulting in 


Cov{fa,Tb) _ 1 _ 1 

V^Var(^)JVax(^ A -|-1 N 


(59) 


Note that the joint probability distribution for 
two partial time delays was obtained in Ref. |160) 
(Eq. 23). 

A systematic expression of the moments (t^) of 
the Wigner time delay for /3 = 2 was given more 
recently by Novaes m\- We see that the vari¬ 
ance diverges for A/3 ^ 2, and the third cumu- 
lant for A/3 ^ 4. Writing the denominator of the 


fourth cumulant as A‘*(A-|-3)(A-|-2)(A-|- 1)^(A — 
2//3)^(A — 4//3)(A — 6//3), we see that (r^)c di¬ 
verges for A/3 ^ 6. These observations suggest the 


power law tail ^-'i-PN/ 2 ^ conjectured in 

Ref. [HI] for /3 = 2 on the basis of some heuristic 
argument involving resonances [i.e. identifying the 


tail of with the one of Eq. (49)]. 

Using a Egdeworth expansion, Mezzadri and Simm 
concluded that the distribution weakly converges 
towards a Gaussian form as A grows. The full 
distribution of the Wigner time delay for A 1 
was however shown to present a richer structure in 
Ref. |182j where the large deviations were studied 
in detail, leading to the behaviours 




for r ^ 0 (60) 

_ g-(/3ivV8)(iVx-i)= for 1/A (61) 
~(At-1)-2-W2 

for At - 1 > y/(2/A) In A (62) 


sketched on Fig. [^ In particular, the transition 
between the sharp Gaussian peak and the power law 


tail was shown to be related to a phase transition 
in underlying the Goulomb gas m- 


Gaussian 



Figure 4: Wigner time delay distribution in the limit N 
1 ; sketch of the rescaled distribution PnIs) = Td = 

ST^), where = h/{NA) is the dwell time. The small 
curves are sketchs of the optimal distributions of eigenval¬ 
ues of Q~^ with the constraint that Tr{Q} is fixed. The 
transition between the sharp Gaussian peak and the power 
law tail is associated with a phase transition in the density 
of eigenvalues. Figure from JTMf . 


Few remarks : 


• A systematic analysis of the “moments” 
(Tr{Q"}) within RMT was carried out by sev¬ 
eral authors : |133111341 l60l 1141] . 


The determination 
and the covariance 
also m) has involved 
tween the two 


of the variance dll^ 
(116) in Ref. [HI] (and 
the correlation be- 
linear statistics Tr{Q} and 


Tr { }. A systematic analysis of the covari¬ 

ance Gov( Tr {Q”} , Tr { Q"*}) was carried out 
by Gunden recently [60], based on the recent 
work |59j where a general formula for the co- 
variance of two linear statistics was derived 
with the Goulomb gas method. More recently, 
the authors were also able to further anal¬ 
yse the correlations (Tr {Q"i} Tr {Q”^} • • • )c 
in Ref. ETl. 


4.1.2. Energy correlations 

The knowledge of the correlation function 

(Tw(e)Tw(e'))c = (T-w(e)Tw(e')) - ('rw(e)) (Tw(e')) 
is of importance and has practical applications : 
for example, this information was used by Polian- 
ski and Biittiker in order to study the effect of ther¬ 
mal fluctuations on the non-linear conductance of 
a quantum dot m- The correlation function was 
obtained by Lehmann et al. m by a random ma¬ 
trix analysis in the orthogonal case (/3 = 1), within 
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the HA of random matrices (see also 


(rw(£)Tw(£')), 

i'Twf 


4 1 - juJTd)'^ 

[1 + (ujTdrf 


(63) 


where ui = e — e'. For e = e' we recover the cu- 
mulant (53) at leading order 1/A^^. Eq. ( [6^ shows 
that correlations occur on the scale l/r^ ^ NA 
much larger than the mean resonance spacing. 0 
We have introduced 

27r 

Td = (tw) = ^ = m/N , (64) 


the dwell time for an electron in the cavity (Eq. 
is valid for strongly overlapping resonances, 1 /t^ ^ 
A, i.e. A ^ 1). The result (63) was later gener¬ 
alised [55] in order to include parametric correla¬ 
tions and describe the crossover between orhtogonal 
and unitary cases. 

Time delay correlations can also be determined 
within the stochastic approach : the model de¬ 
scribed in Refs. [U 11481 l3T] and used by Polianski 
and Biittiker [1491 115011151] provides some infor¬ 
mation about the energy and magnetic field depen¬ 
dence of the 5-matrix correlator, which reads [TM] 

{Sabie, B)S*^{e\ B')) = Sac^bd ^e-e' + Sad^bc ■ 

(65) 

and are the (zero-dimensional) analogues of 
the Diffuson and the Cooperon appearing in the 
diagrammatic approach for weakly disordered met¬ 
als [1] , [^ given by 




NTd l/r®,<g’ - iuj ' 


The two characteristic times 


1 1 1 

- —-1- 


( 66 ) 


(67) 


In a ballistic QD of size L, the Thouless energy, which 
sets the energy scale for correlations in a closed QD, is 
= vpjL, where vp is the Fermi energy. The open QD 
is characterised by the number of channels N = k^w/Tr, 
controlled by the width w of the contacts. In this case, as 
we have seen, the relevant scale is a different Thouless en¬ 
ergy [13 We have e^^^/sTh '^/L < 1. 

This can be understood from the Fisher and Lee relation 
m between scattering matrix and Green’s function. 

Note that ( |65| l only satisfies unitarity at leading order 
in the orthogonal case : one gets {{S^S)aa — 1) = = 1/N 

instead of zero. This can be corrected by using the expression 

HI {Sab{£)S*^{E')) = [&aa&bd + (2//3 - l)(5ad<56c] with 
[Ar + 2/d - 1 - 2i7ra;/A]“b 


combines a contribution describing the escape rate 
l/Td from the cavity and dephasing due to the 
magnetic field. The magnetic dephasing rate is 
l/'^B = (^iF^/Surf)(<i)/(/)o)^ where vp is the Fermi 
velocity, $ = Surf B the magnetic flux through the 
cavity and (j)Q = h/e the quantum flux, i is the size 
of the cavity in the ballistic case or the elastic mean 
free path in the weakly disordered (diffusive) case. 


As a simple application of the formula we get 
(tw(£)) — Td{N‘^^Q + N^q ), where we have used 
-\de^e-e' = fh ■ We deduce (tw(£))/ r^ ~ 

1 -I- 0{1/N), whose leading order term coincides 
with ( |5^ . 


at 

ex- 


The analysis of the correlations requires 
additional information : simply apply¬ 

ing Wick’s theorem with the correlator 
(j^, we get (rw(£,H)rw(£',^'))!i°''''"’'‘“^ - 

-I- |^E_e/|^), which leads, 

zero magnetic field, to the incomplete 
pression (tw(£, 0)tw(£', ('^w)^ - 

{2/N^) [l -I- {iOTd)^] ^. This calculation, known 
as the “diagonal approximation” in the context of 
semiclassical methods, disagrees with ( |6^ : not 
only the energy dependence differs but also the 
value at a; = 0 is half the correct result [compare 
with the variance (53)]. This emphasizes the 
importance of non Gaussian contributions to the 
correlator, which are taken into account in the 
more precise expression of the four-point correla¬ 
tion function {SabS*^SijS, 
of Ref. |3T| (see also |57]1. 


j,;) given in Appendix B 


This question has been much discussed in semi¬ 
classical approach : this discrepancy was under¬ 
lined by Lewenkopf and Vallejos [120] . Kuipers and 
Sieber m identified the nature of the contribu¬ 
tions (“trajectory quadruplets”) to be added to the 
diagonal approximation in order to recover the ran¬ 
dom matrix result (63) (doing so they have recon¬ 
ciled the two semiclassical approaches for a Wigner 
time delay analysis : the periodic orbit expansion 
based on the relation with the density of states 
and the scattering approach involving trajectories 
entering and leaving the system). The expression of 
the correlator, determined earlier within a semiclas¬ 
sical approach by Vallejos et al. [189] (for B = B'), 
is the sum of two contributions which can be identi¬ 
fied as Diffuson and Cooperon contributions. Rein¬ 
troducing the effect of the field mismatch, we get 
the structure obtained within the HA of random 
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matrices in 

Ref. mUl] : 


the chiral-unitary class (AIII) in the N = 1 channel 
case : 

(rw(e, 

S)rw(e',S'))c 

(tw)^ 

(68) 

.^32™(r) = r-2e-i/" (69) 

in appropriate units. This should be compared with 

2 

/ + 

1 

the expression (41) for the unitary case (class A) : 


l[i/4 + "4 

[l/T2+a;2]^j 

We see that the exponent of 


where we recall that th = 27r/A. This expression 


now agrees with (631 for N. 

Few remarks : 

• The time delay may be represented in terms of 
as tw(£:) = (l/Af)Eara/[(£ - ea)2 +r2/4], 
where the sum runs over the resonances Ena 
IM] . This establishes a connection between 
Wigner time delay fluctuations and Ericson 
fluctuations of the cross-section for strongly 
overlapping resonances T > A iMiina- 


the tail is increased from 2 to 3 when the chiral 
symmetry is broken, as for the ID disordered case 


reviewed above : compare (37) in the presence of 


the chiral symmetry and (25) when the chiral sym¬ 
metry is broken. 


More recently, the result (42) for the joint distri¬ 


bution of proper time delays was extended to the 
other symmetry classes in Refs. [1241 I162j . This 
question is reviewed in another article of this spe¬ 
cial issue |125j . 


• Semiclassical methods were widely used in 

order to analyse Wigner time delay correla¬ 
tions ; correspondence 

between the random matrix and semiclassical 
approaches holds for strongly overlapping res- 
onaces I/t^ I/th ~ A (i.e. N :$> 1) |189j . 

• More recently, an improved semiclassical ap¬ 
proach for the analysis of the Wigner time 
delay statistics was developped by Kuipers, 
Savin and Sieber m, who carried out a di¬ 
agrammatic calculation of the moments of the 
Wigner time delay. A semiclassical derivation 
of the moments of the time delay was also pro¬ 
posed by Novaes |142j . 

4-2. Other symmetry classes 

As we have already mentioned, the three Wigner- 
Dyson symmetry classes (denoted AI, A and All 
in the Altland-Zirnbauer classification) were com¬ 
pleted by three chiral classes (chiral orthogonal 
BDI, chiral unitary AIII and chiral symplectic CII) 
and four Bogoliubov-de Gennes classes (C, Cl, D, 
Dill) [wiiHirfu. The new symmetry classes have 
attracted a lot of attention during the last years 
in relation with topological insulators [Ml I163j and 
topological superconductors [mds] (see also [158] 
were a classification of the various types of topolog¬ 
ical insulators was provided for the different sym¬ 
metry classes and dimensions). 

The analysis of the Wigner time delay distribu¬ 
tion for different symmetry classes has been initi¬ 
ated by the work of Fyodorov and Ossipov |5T] for 


5. Generalised concepts : injectance, emit- 
tance, etc. 


The concept of scattering matrix is central in the 
Landauer-Biittiker approach. For this reason, the 


Krein-Friedel relation (16) plays a very important 


role as it allows to express the density of states 
(DoS) of the open conductor in terms of scattering 
properties. In the following, we will simply write 


-(£) = A:Tr{C} 


2i7r 


de 


2 nH 


(70) 


where we have neglected the term Tr — 5^} /e 
in Eq. (16) (this is justified in the metallic regime). 
Such relation allows to characterize the amount of 
charge injected in an open coherent conductor, a 
crucial tool which has been used by Biittiker in or¬ 
der to describe screening properties. The out-of¬ 
equilibrium situation where the conductor is con¬ 
nected to several contacts (terminals) with differ¬ 
ent chemical potentials however requires to iden¬ 
tify the contributions of each terminal to the DoS, 
which has led Biittiker to introduce several gener¬ 


alisations of (70) as the partial DoS, the injectance 


and the emittance. These concepts have been used 
to develop a theory of non-linear transport |38l EH] 
and AC transport [38l Ell Ell EH 11621 in coher¬ 
ent conductors (see the reviews Ha 13111111115]). 


5.1. Partial DoS, injectance and emittance 

Let us consider a multi-terminal structure, whose 
scattering properties are characterised by a basis of 
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stationary scattering states '!/'e,a(a^) where the in¬ 
dex ot labels the terminals. We assume that the 
terminals support each a single conducting channel 
(i.e. contact wires are effectively one-dimensional), 
what simplify the analysis (the generalisation to 
contacts with several channels is straightforward). 
Furthermore, this would allow to illustrate the dis¬ 
cussion by explicit formulae by considering the case 
of metric graphs for which explicit construction of 
the scattering matrix is possible [18411178j . 

Biittiker introduced the concept of '’'’partial den¬ 
sity of states’' [SHI 09] 



measuring the contribution to the DoS (70) of par¬ 


ticles incoming from terminal /3 and outgoing at 
contact a. Another quantity which appears in 
Biittiker’s work is the "injectance ", obtained by 
summation of 0 over the first index 




(72) 


It provides the contribution to the DoS of the scat¬ 
tering states incoming from terminal a. Summation 
over the second index leads to the " emittance” 

ka(e) = ^^a/3(£)-^ (^5^) (73) 

^ \ / aCK 


which corresponds to the contribution of particles 
outgoing at terminal f3. Injectance and emittance 
are related by magnetic field reversal. 


• These concepts have been further generalised 
in order to deal with local properties, leading 
to the concept of partial local DoS, injectivity 
and emissivity jSSlIEHj • This can be understood 
from the relation 


where (x) is the potential (see Appendix of 
Ref. [10] or Refs. |175[ 11781 IlSOj b This leads 
to introduce the injectivity 


Va{x’,e) 


1 

2i7r 




5 S 

5^ (x) 


ack 


IV'E.a(a()|^ 

(77) 


which measures the contribution of the scat¬ 
tering state to the local DoS y(x-,e) = 

J2a^oi(x'j Similarly, Biittiker intro¬ 

duced the concepts of partial local DoS and 
emissivity [551155] . 


• The Fisher and Lee relation m = 

—Sa 0 + iy/VaVp G^{a,/3;e), where Va is the 
group velocity in terminal a and where 
the retarded Green’s function is measured 
at the two terminals, makes the functional 
derivation transparent : dSapfS^(x) = 
iy^FlTFa G^(a, x; e)G^(a7,/3; e). Cf. also the 
discussion in [178] . 


• An illustration : the knowledge of the injectiv¬ 
ities allows to express the density of electrons 
in the conductor out-of-equilibrium 


n(x)=^ J de f{e - eVa)iya{x’,e) 


(78) 


= Va{e]-B) , (74) 

which follows from S{—B) = S{B)'^. Obviously 
Hs) = = '^Val5i£) ■ (75) 

O' a 

Injectance and emittance allow for a preselection 
and a postselection, respectively, when identifying 
the contribution to the density inside the conductor 
from particles passing through it [55] . 

Few remarks : 

• All formulae can be straightforwardly gener¬ 
alised to the case of multichannel contacts by 
adding some traces over channels. 


where /(e) is the Fermi function and Va the 
electric potential at contact a. 


• The introduction of the injectivities allows to 
avoid the approximation made in (70) by ne¬ 
glecting the term Tr{5 —5'^} /e of Eq. (16), 
and provides a better definition for the partial 
DoS, injectance and emittance. Introducing a 
uniform potential (x) = U inside the con¬ 
ductor (but not in the contacts), we can write 


(79) 

where /q^ • • • denotes integration inside the 
conductor [the boundaries of integration are 
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given by the place where the scattering states 
are matched with plane waves in order to de¬ 
fine scattering amplitudes, like in Eq. ([^]. In 
other terms [17811180| 


Therefore, instead of (721 which involves a high 


energy approximation, a more rigorous defini¬ 
tion of the injectance should be 


.(e) = - 


1 

2i7r 




95 


(81) 


The remark also holds for partial DoS and 
emittances, Eqs. ( 71|73 |. The DoS takes the 
form i/{e) = — (2i7r)“"9£/ lndet5. 


5.2. A generalisation of the Feynman-Hellmann 
theorem for a continuous spectrum 
Biittiker’s idea of relating the spectral properties 
to the scattering matrix is not limited to the DoS 
and the local DoS and can be generalised. This may 
be viewed as an extension of the famous Feynman- 
Hellmann theorem which applies to bounded prob¬ 
lems with discrete energy spectra. Let us first re¬ 
call this well-known theorem : consider a Hamil¬ 
tonian H characterised by the discrete spectrum 
{sn, I'*/’«)}■ The determination of the diagonal ma¬ 
trix elements of an observable X = —dH/df, where 
/ is a conjugate force, does not require the knowl¬ 
edge of the eigenvectors but only of the eigenval¬ 
ues : {'ijjn\X\'ipn) = —dsnldf. In Ref. |180j . we 
proposed that an extension of this theorem to the 
case of Hamiltonians with continuous spectra, char¬ 
acterised by scattering stationary states {IV’e.a}}) 
is 

^ ^ . (82) 

This relation was proved in |178[I180] for various ob¬ 
servables in metric graphs : for example, choos¬ 

ing the local density X —)• p{x) = |a;)(a;| with 
/ —>■ we recover ( [76| . In |18()j . we con¬ 

sidered also the case of the current density la in 
a wire (a) of the graph, which involves a deriva¬ 
tion with respect to the magnetic flux (fa along the 
wire {ipe,a\la\'ifs, 0 ) = { 2 m)-^{S^dS/d(l)a)^^- A 


Note that the expression ( |82| l does not account for the 
contributions of BICs if present, as the DoS [Tsoi. 


trace over indices leads to the well-known formula 
obtained by Akkermans, Auerbach, Avron and 
Shapiro [3] : (7)^ = (2i7r)“^9,^ lndet5. Note that 
similar considerations were used in order to analyse 
current fluctuations at equilibrium in Ref. |172j (a 
formula for current correlations out-of-equilibrium 
was obtained in |I8()j l. 


5.3. Statistical analysis in quantum dots 

In the case of chaotic quantum dots with several 
contacts with perfect couplings (Fig. [^, Brouwer 
and Biittiker [28] obtained the first two moments 
of the partial DoS within the ASA of random 
matrices : EH 


1 

\NpN, 

1 


[NpN, 

'll 

A 

N2 

N 


N J 


(83) 


where is the number of open channels at contact 
p and N = We deduce (i7„) = {p_a) = 

(Na/N) as expected in an ergodic device. 

In the unitary case (/3 = 2), the covariances 
are |55| 


Cov(j 


fllA 5 ‘^p(T 


) = 


[NpN, 

fs 


iV2 


+ A j 


3 



(N,Na 


2 

{ A2 

A 

A2 

A \ 


(84) 


The correlator for the orthogonal case (/3 = 1) is 
obtained by adding similar terms obtained by per¬ 
mutation p -f-)- (7 in order to fulfill the symmetry 
Vpa = Vap (for H = 0) ; the symplectic case receives 
additionally a factor 1/4. Eq. (84) shows that meso¬ 
scopic fluctuations are of order Svap {vap)/N ^ 
1/(7VA). 

The correlations of injectances and emittances 
are also of interest. We deduce (for /? = 1, 2) 




1 

(AA)2 


NpN, 

iV2 



/2 A 2NpN,' 

U ) iV2 

(85) 


Brouwer and Biittiker considered the dimensionless AC 
conductances = Ga.^{0) — ia; + O(aj^), which, 

in the unscreened limit, can be related to the partial DoS 
: limc->oo = 27r/d£(-as/) J.a;3(e). 
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and 


Cov(i 


fj-’1^1/ 


) = 


{NAy 

2 N^N, / 2 

7V2 + ' " 


N2 n 


( 86 ) 


The interchange of the two contributions between 


(851 and (861 is analogous to the exchange between 


Diffuson and Cooperon in weakly disordered metals 
|183j . These expressions allow us to quantify the 
difference between injectance and emittance 


((^a-Tj^)= 2- 


/3 


2 iV„ 
{NAy N 


1 - 


N 


(87) 

The difference vanishes in the orthogonal case {/3 = 


1), as a consequence of the symmetry (74). This 
quantity has found some application in |159j . 


6. Non-linear transport in coherent conduc¬ 
tors 

This section gives a brief presentation of the the¬ 
ory of non-linear transport in coherent conductors 
proposed by Biittiker and Christen [38l[5llll3], il¬ 
lustrating the interest of the quantities introduced 
in the previous section. We consider the case of 
conductors of mesoscopic dimensions such that they 
can be considered in the ergodic regime (like quan¬ 
tum dots), which slightly simplifies the presenta¬ 
tion, although the theory was developed in a more 
general context in Ref. [3SIISI]- 



Figure 5: A conductor of mesoscopic dimension connected 
to M = 3 macroscopic reservoirs and capacitively coupled to 
a gate. When the potential at contact a is raised, current is 
injected in the conductor from this contact (blue arrows). 


6.1. Non-linear conductances 

We consider a multiterminal mesoscopic struc- 
ture (Fig. [^. For small voltages Va —t 0, current 


in terminal a can be written under the form of an 
expansion 


c = 



gafi Vfl + 
/3 



^ ^ 5a/37 VaFy -l- • • • 
P 


( 88 ) 

where 2^ is the spin degeneracy. gajS is the dimen¬ 
sionless Landauer-Biittiker linear conductance m 


9ct^ — 


de 



(89) 


f{e) is the Fermi distribution and gapis) = 
— Tr{5(^^(e)iSQ,/3(£)} the zero temperature 
dimensionless conductance, Na being the number 
of conducting channels in contact a. 

The Landauer-Biittiker formula = 

{ 2 se/h)J 2 i 3 J —eVp) gap{e) already produces 
a contribution to the non-linear conductance 

glp-y = j • (90) 

Biittiker emphasized the importance of screening 
effect which produces a second contribution that 
we now determine. The first step is to analyse the 
electrostatic inside the conductor. 


6.2. Characteristic potentials 


Let us recall the main ingredients involved in 
the description of screening in good metals |196] . 
which Biittiker adapted to deal with the non¬ 
equilibrium situation. These analysis is based on 
three ingredients. (*) When charges are intro¬ 
duced in the conductor, the total density of carri¬ 
ers Sntot{x) = Suextix) + Suindix) is the sum of the 
density Snext(x) of charges injected in the conduc¬ 
tor and the density Snind(x) induced by screening. 
The concept of injectance introduced above allows 
us to write the number of injected charges when 
the conductor is out-of-equilibrium. An expansion 
of Sne^t(x) = n(x) — nequii(x) where n(x) is given by 
(78) and nequii(a^) is the equilibrium density, gives 


M 


^^ext — 


'QD 


da;i5next(a;) — eVai^aisr), (91) 


a—1 


where the summation runs over the M contacts. We 
assume zero temperature for simplicity ; the finite 
temperature formula involves an additional convo¬ 
lution with a Fermi function J de ( — df /de) ■ ■ ■. 
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{ii) The induced density is related to the potential 
energy ‘W (x) by linear response theory 

6 nind{x) = — J dx'n(x, x') ^(x') (92) 

where n(a;,a;') is the (static) density-density corre¬ 
lation function of the non interacting electron gas 
(Lindhard function [ddj 1. Assuming ergodic prop¬ 
erties, the potential energy can be considered uni¬ 
form in the conductor (x) ~ U. Thomas-Fermi 
approximation relates the response to the DoS : 


Suind = / da;(5nind(a::) ^ -v^sf) U . 

JQD 


(93) 


(Hi) The potential energy and the density are 
related by the Poisson equation (x) = 

—dne'^ Sritotix) which encodes the Coulomb inter¬ 
action in the conductor. The fact that (92) involves 


the potential (x) and not the electrostatic poten¬ 
tial related to the external density Sriextix) makes 
the approach self-consistent. For a uniform poten¬ 
tial, the Poisson equation is replaced by 

C {U - eVo) = e^{6next + (Juind) (94) 

where C is the capacitance of the conductor and Vq 
the electrostatic potential of the gate, labelled with 
index 0 (Fig.[^. Excess charge inside the conductor 
occurs when the system is driven out-of-equilibrium 
by the external potentials. At linear order, it is 
convenient to decompose the potential over contri¬ 
butions of the different external voltages 


M 

U ~ Ua eVa 

a=0 


(95) 


where the ^^characteristic potentiaF Ua measures 
the response of the potential C7 to a shift of the 
external voltage Va [381E] (one should not forget 
the response of the gate voltage controlled by Ug). 
Injecting ( 91|93|95' ) in (94) provides the character¬ 
istic potentials 


C'm ^a(^F) 
C iy{sF) 


for a G { 1 ,..., M} and 




(96) 


(97) 


We have introduced the '’^mesoscopic capacitance 
1 1 1 


C„, 


c a 


(98) 


combining the "geometric capacitance''^ C and the 
"quantum capacitance" 

Cq = e^ vIef) ■ (99) 

We can check the sum rule 

M 

^Uc = l. (100) 

a=0 


The ratio yint = Cp^jC in Eq. (96) may be viewed 


as an interaction constant as it measures the effi¬ 
ciency of screening (yint — 1 for efficient screening 
in a good metal and yint ^ 1 for weak screening). 
The response of the gate ug = 1 —7int is thus impor¬ 
tant when screening is weak, however the presence 
of the gate (or the surrounding medium) can be 
forgotten in good metals. 

6.3. Non-linear conductances 
We can now determine the non-linear conduc¬ 
tances as follows. The redefinition of the elec¬ 
trostatic potential inside the conductor leads to a 
modification of its scattering properties which de¬ 
pends on the potential, what contributes to the 
non-linear response : gapisp) —t gapi^F — U) ce 
gapi^p) - g'api^p) ^7 6 ^ 7 . After symmetrisa- 

tion with respect to indices /Joy and reintroduc¬ 
ing the Fermi function, we get the expression of the 
non-linear conductance [HI I149j 




( 101 ) 


X [9'o,p(£)^P'r-g'ap{£)ui-g'a^{s)up] 1 
where the first term comes from the noninteract¬ 


ing theory, Eq. (90). For consistency, the T = 0 
expression ([%]) should be replaced by 


Jds{-def) Va{e) 


( 102 ) 


C/e^ + J de (-4/) v{e) ' 

The second term of the denominator is the h- 
nite temperature quantum capacitance Cq = 
e^ J ds{-dJ)F{e). 

Few remarks : 

• As emphasized above, all quantities involved in 
Eqs. ( |101|102| ), g OLp{e), v{e) and Vaie), are ex¬ 
pressed in terms of the iS-matrix. The deriva¬ 
tive of the conductance g'apis) can also be re¬ 
lated to the sensitivies introduced by Gaspar- 
ian, Christen and Biittiker [SHI : 




dSap 
de 


dS. 


0.(3 


de 


S, 


0.(3 


(103) 
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(note that Ref. |HH1 rather introduced a local 
version of the sensitivies). 


and specifically the asymmetry in magnetic 
field pT^rTMITMl I^I^. 


A statistical analysis of the non-linear conduc¬ 
tance (1011 in chaotic quantum dots requires 
the statistical properties of the conductance’s 
derivative provided in Ref. [32] and those of 
the characteristic potentials, i.e. of the injec¬ 
tivities, obtained in [21], see Eqs. (83|84). 


6.4- Recent developments 

• As mentioned, the Biittiker’s scattering for¬ 
malism can be extended to discuss conductors 
which are not in the ergodic regime [381 IhTl 
133] . The equivalence between this formalism 
and the non equilibrium Green’s function ap¬ 
proach for non-linear transport was established 
by Hernandez and Lewenkopf |l()2j (see also 
[183]). 

• An interesting aspect was identified by 

Sanchez and Biittiker [159] and Spivak 
and Zyuzin |169j : the interacting part 
of the non-linear conductance = 

“(l/2)[5a/3(^F)% + g'a-y{^F)uii\ is not Con¬ 
strained to any specific symmetry under re¬ 
versal of the magnetic field, contrary to the 
non-interacting part = {l/2)g'^p{£F)5p^. 
This is due to the asymmetry of the char¬ 
acteristic potentials Ua (i.e. of the in- 
jectance Va)- The analysis of the asymme¬ 
try of the non-linear conductance under mag¬ 
netic field reversal [giii(H) — giii(—H)]/2 = 
— {Ci_i/C) giiisp) {fi — Ti)/(2 j^) was thus pro¬ 
posed as a new way for probing electronic in¬ 
teractions in mesoscopic structures (the pa¬ 
rameter 7 i„t = CfijC introduced above was 
measured in [9] for a mesoscopic ring). Spi¬ 
vak and Zyuzin borrowed arguments from di¬ 
agrammatic techniques in order to analyse the 
low magnetic field regime |169j . Sanchez and 
Biittiker |159j proposed a random matrix ap¬ 
proach in the unitary case (strong magnetic 
field), which was later improved by Biittiker 
and Polianski who extended the statistical 
analysis in order to describe the crossover from 
orthogonal to unitary cases and include ther¬ 
mal effects |14911150] (for more details, cf. the 
excellent review article [151j l. Note also the 
study for a ring made of strictly ID wires m- 

• Several experimental groups have analysed the 
non-linear transport in mesoscopic structures 


• The case of diffusive wires is considered in 

Ref. [T33] . 


7. AC transport 

The search for fast control and manipulation of 
charge in coherent conductor has stimulated many 
developments (see the reviews [86l[23]). Biittiker, 
Pretre and Thomas have proposed a theory of 
time-dependent response in coherent conductors 
gTiiiHiiin] based on the scattering approach and 
a Hartree-Fock treatment (with Thomas-Fermi ap¬ 
proximation) of electronic interactions. In order 
to emphasize few ideas, we will restrict ourselves 
here to the case of the “quantum RC circuit’ which 
was studied experimentally in the integer quantum 
Hall regime during the PhD of Gabelli [H3 [T3| 157] 
(Fig. §. The RC circuit was later studied from 
the perspective of an emitter of electronic wave 
packets in the PhD of Feve [73] and as such is a 
building block for “electron optics experiment” (see 
the review [23] and article m of this special is¬ 
sue). These latter developments have mostly con¬ 
sidered the RC circuit in the integer quantum Hall 
regime, such that the current is carried by a single 
edge state along the boundary of the system, when 
one conducting channel is open at the constriction 
(quantum point contact, QPC). 

Below, we do not consider the integer quantum 
Hall regime but we rather discuss the situation 
where the RC circuit of Fig. [^ is submitted to a 
weak magnetic field, such that the electron dynam¬ 
ics inside the cavity is chaotic. This justifies a ran¬ 
dom matrix approach. In this case, the opening of 
the constriction controls the number of channels N. 



Figure 6: The quantum RC circuit : a coherent conductor 
patterned in a 2DEG (blue) is closed by a quantum point 
contact (QPC) controlled by side gates. A capacitive cou¬ 
pling to a top gate allow to analyse its AC response. Figure 
from Ref. J91f . 
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7.1. The AC response 

By using arguments similar to the ones described 
in Section Biittiker and coworkers have shown 
that the Coulomb interaction can be taken into ac¬ 
count by adding the impedance of the non interact¬ 
ing electrons, 1/Go{uj), and the classical impedance 
of the capacitance Z{uj) = 1/G'(w) = l/G'o(<j-') + 
l/(—iwC). The admittance is [TT] 


G'o(u;) = J de Tv {l - S\e)S{e + oj)} 
f{s) - /(e + w) 


(104) 


where 2s denotes the spin degeneracy. One can 
write the low frequency expansion of the impedance 
under the form given by the elementary laws of elec¬ 
trokinetics 


Z{oj) 


+ R^ + 0{u;), 


(105) 


where the ^‘‘mesoscopic capacitance’’’’ and the 
“charge relaxation resistance” Rq carry some infor¬ 
mation about the dynamics of electrons in a quan¬ 
tum coherent regime, as they control the relaxation 
time Tuc = RqC^. The presence of the mesoscopic 
capacitance (98) in the AC response comes from 
the a; —)■ 0 expansion of Go(w), whose first terms 


involve the Wigner-Smith matrix (12). Hence, this 


latter controls both the quantum capacitance (99) 




2se2 


Tr{Q} 


(106) 


and the charge relaxation resistance 
h Tr{g^} 

2sX2e2(Tr{Q})2 


(107) 


In these equations, the Wigner-Smith matrix is 
taken at Fermi energy ep. 


7.2. The charge relaxation resistance 
In this section we forget about spin degeneracy 
factor 2s and base the discussion on the expression 


, h Tr{g^} h Y^grl 
2eMTr{Q})2 2e^ ' 


(108) 


A first important observation is that the resistance 
belongs to the interval : 


h 

2N^ 



(109) 


In particular the A^ = 1 channel case leads to the 
universal value Rq = h/ (2e^) [48] (this is a property 
of a coherent conductor, cf. remarks closing the 
section). 


Case N = 2 —. Assuming a random matrix de¬ 
scription of the quantum dot, we can deduce the 


distribution of Rq starting from (42). The calcula¬ 


tion only involves a double integral 


P2 


{rq)=A^{2rq-l)(P-^^/^{l-rqf ( 110 ) 


where rq = {Ne^/h)Rq is the dimensionless 
charge relaxation resistance per channel ; Ap = 
2^+YB{f3 + l,{f] + l)/2) is a normalisation con¬ 
stant. We deduce the mean value 


{Rq) = 


h 


( 111 ) 


which is surprisingly independent on the symmetry 
class. 


Remark : The distribution of Rq for A^ = 2 
channels was analysed by Pedersen, van Langen 
and Biittiker in |146j . however these authors in¬ 
troduced a slightly different averaging procedure : 
following [32] the authors weighted the joint dis¬ 
tribution of proper times by the DoS, correspond¬ 
ing to a “canonical” averaging. They deduced 
_ i)(/3-i)/2(x — rq)^~^ where 
the normalisation is Bp = 2^/i?(/3, (/3 -|- l)/2). As 
a result the mean resistance depends on the sym¬ 
metry index = (3/8) h/e^ for /3 = 1 and 

= (5/14) h/e^ for /3 = 2. 

Mean value —. Biittiker, Pretre and Thomas [48] 
already pointed out that the charge relaxation re¬ 
sistance should scale with the number of channels 
as Rq ~ 1/A^, which corresponds to the addition 
of resistances in parallel, although the precise pref¬ 
actor was not given. The fact that the one chan¬ 
nel case (spin polarised) leads to the universal half 
quantum of resistance Rq = hl{2e^) has produced 
a certain confusion for the large N case, between 
Rq = h/(2A^e^) (incorrect) and Rq = h/{Ne^) (cor¬ 
rect). 

The hrst precise calculation of the mean value 
{Rq) in the A^ 1 limit is contained in the study 
of the complex admittance G(w) by Brouwer and 
Biittiker [28] within the ASA of random matrices : 
the low frequency expansion of (G(a;)) (Eq. 12 of 
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Ref. 


leads to 


(this coincides with the leading order term of cumu- 


{R,) = 






-+0(iV-2) 


( 112 ) 

(without spin degeneracy 2s). The second term can 
be interpreted as a weak localisation correction. 

A first review paper by Biittiker |39j incorrectly 
quoted this result by introducing a spurious factor 
of 1/2 (see Eq. 17 of Ref. [39]). In a second review 
PP] . Biittiker gave the correct result [Eq. 9 of this 
reference] 

(«,) = ^ (113) 

for a non perfect contact with NT ^ \ (without 
spin degeneracy), where T is the transmission prob¬ 
ability through the contact. Biittiker refered to a 
private communication with Carlo Beenakker and 
some unpublished work with myself that I 

describe in Appendix A This value coincides with 
the DC resistance of the constriction, what was ex¬ 
pected in the large N limit. 

We can summarize this discussion in the follow¬ 
ing table giving the dimensionless charge relaxation 
resistance per channel Vq = {Ne^ /h)Rq : 

{rq) = {Neyh){Rq) 


N = 1 
N = 2 
A^» 1 


1/2 

2/3 

l + {2/P-l)N-^ 


Fluctuations —. As the quantum capacitance is di¬ 
rectly proportional to the Wigner time delay, its 
statistical properties were determined in |182) for 
a contact with many channels. More recently, we 
have studied with Grabsch m the statistical prop¬ 
erties of the charge relaxation resistance, which in¬ 
volves the analysis of the ratio of two “linear statis¬ 
tics” X)a'^a/(So'^a) ' ^e have obtained the mean 
values (Cq) ~ e^/A and (Rq) ~ h/{Ne^) (without 
spin degeneracy), in correspondence with the DC 
resistance of the QPC, as it should. The variance 
of the capacitance is given by |1171 |2H1 118211135L jPD] 


{5Cl) 4 

{Cqf 


(114) 


Although the geometries are different, we can compare 
this expansion with those of the DC resistance of a coherent 
quantum dot closed by two constrictions with and N 2 
channels [13] : {Rdc> = Kdass [l -t (2//3 - + ■■•] 

where i?c:lass = (1/Afi -|- I/IV 2 ) and N = + N 2 . 


lant (531 for TV ^ 1). The variance of the resistance 

(K) 


{R, 


8 


(115) 


and the correlations with the quantum capacitance 

(SCqSRq) 


^{5Cl){5Rl] 


— + 1 / 


(116) 


were found in Ref. m, where we have also analysed 
the large deviations. 

7.3. Final remarks 

• The multiterminal version of the general for¬ 
mulae reviewed in the section was provided in 
Refs. 1^51 HP] . 

• The finite temperature formulae are obtained 
by introducing convolutions with Fermi func¬ 
tions J de ( — df/de) : the traces over channel 


indices in Eqs. (1061107) are then affected by 
additional integrations over energy. 

In particular, using the representation Cq = 
(e^/h) f de {—def) Tr {Q(e)} with the correla¬ 
tor ( 681 , the model leads to the following ex¬ 


pression 

{SCq{B)6Cq{B')) 


(O 

1 


~ J dujST{aj) ( 117 ) 


In [(I -F (uTsif) (1 -I- (ujTc^f)] 


where St(uj) is a (normalised) thermal func¬ 
tion of width T which arises from the con¬ 
volution of the two Fermi function’s deriva¬ 
tives. 1^ For Trdi^ 1 this leads to the decay 
{SCqiBf) / {Cqf ~ [K/(^iV2)] (TTd)-2 where 
K = — dxF'{x)/x ~ 0.365 (the function 
F{x) is defined in the footnote). The depen¬ 
dence in the Dyson index jS can be replaced by 
the precise magnetic field dependence straight¬ 
forwardly. 

Note that in the metallic limit (e^/C A, 
i.e. C ^ Cq) the fluctuation of the mesoscopic 
capacitance is ~ (C/ {Cq})'^ SCq. 


23 StO = F{uj/{2T))/{2T) with F{x) = (a;cotha; - 
1)/sinh2 X. 
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• A more precise discussion of the capacitive re¬ 
sponse of quantum dots has been provided in 
Ref. [5], beyond mean field and also analysing 
the limit of a weakly transmitting constriction 
in the presence of Coulomb blockade. 

• Nigg and Biittiker have considered in Ref. [139] 
the case of a single spin polarised channel. 
They showed that strong enough dephasing in 
the cavity can induce the transition from the 
universal value Rq = h/{2e^) to the DC resis¬ 
tance i?dc = h/{e^T) of the QPC, where T is 
the transmission probability throuhg the QPC. 

• The universal charge relaxation resistance 
Rq = h/{2e^) for a single spin non degener¬ 
ate channel [Ml HU iO] has stimulated sev¬ 
eral works. This value was later obtained 
within more accurate treatments of electronic 
interactions : within a mean-field Hartree- 
Fock analysis |140L I157j improving the original 
Thomas-Fermi treatment, and beyond mean 
field [571 1137] , within Matveev formalism [128] 
of QPC (see also [70]!. These treatments have 
accounted for possible Coulomb blockade ef¬ 
fects. 

The treatment of interaction in a multichannel 
contact was extended more recently in [66] . 

8. Conclusion 

Despite the huge diversity of Biittiker’s work, this 
article has tried to show that a possible Ariane’s 
thread is furnished by the concept of time delay 
and related quantities. Biittiker’s interest in this 
matter goes back to his early fundamental study of 
times in quantum mechanics ; this question has led 
to many possible definitions, from scattering phase 
shifts (Wigner time delay, transmission group de¬ 
lays,...), clock approaches, etc. When randomness 
is introduced, these quantities are characterised by 
non-trivial distributions. I have reviewed in Sec¬ 
tion [^ the statistical analysis for the most simple 
of these fundamental times in the ID situation in 
the presence of disorder. In particular I have em¬ 
phasized the role of symmetry by considering the 
effect of a chiral symmetry. In this case, a new 
observation was the relation between the presence 
of a chiral zero mode and the existence of a limit 
law for the Wigner time delay when the size of the 
disordered region goes to infinity. 

One of Biittiker’s major motivation for Wigner- 
Smith time delay matrix’ analysis comes from its 


central role in the scattering approach for non¬ 
linear (Section]^ and AC (Section]^ coherent elec¬ 
tronic transport. This has led to introduce several 
new concepts (Section which have been anal¬ 
ysed within the frame of random matrix theory : 
Biittiker, with Gopar and Mello, obtained one of 
the first important result in the Wigner time de¬ 
lay statistical analysis [90], motivated by the sta¬ 
tistical analysis of the mesoscopic capacitance of a 
chaotic quantum dot. The first moments of the AC 
response were obtained with Brouwer [28] within 
the random matrix theory and the non-linear re¬ 
sponse with Sanchez and Polianski |159[ 1149] . In 
the present article I have provided a review of the 
known results on time delay statistics within the 
random matrix approach : joint distribution of 
proper times, moments of Wigner time delay and 
its correlation function and statistical properties of 
the partial DoS. 
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Appendix A. Calculation of {Rq) for N ^ 
1 — Mapping to the study of 
absorption in the cavity 


This section reproduces some unpublished con¬ 
tribution due to myself and Biittiker [177] (quoted 
in Ref. [40] but remained unpublished so far). 

Our aim is to estimate the charge relaxation re¬ 
sistance for a one contact conductor in the good 
contact limit NT ^ 1, where N is the num¬ 
ber of conducting channels and T the transmis¬ 
sion probability through the QPC. Our starting 
point is the expression (108) (without spin degen¬ 
eracy) [M]- As we mentioned in the text, the trace 
in the denominator is a measure of the density of 
states (DoS) of the cavity v ~ Tr{Q}//i, there¬ 
fore (Tr{Q}) = 27r/A where A is the mean level 
spacing and the main question left is to estimate 
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Tr { }. We remark that it appears in the expan¬ 

sion of the dimensionless AC conductance (104) 


Goiuj) = N-Ti {S^e - + ^lo)} 


= -iwTr{Q} + — Tr{S2} + ... (A.l) 


Biittiker suggested to use analytic properties of the 
scattering matrix, w —>■ iF, i.e. to map the problem 
of AC transport onto the study of DC transport 
for a cavity in the presence of absorption, whose 
dimensionless conductance reads 


G = Go(ir) = iV - Tr - ^r)5(e + ^F)} 

= FTr{Q}-y Tr{Q2} + ... . (A.2) 


This corresponds to a situation where absorption 
takes place uniformly in space with absorbing rate 
F. In the ergodic regime where RMT holds , this 
also describes an absorbing contact (Fig. A.7).m 



absorbing 

probe 


Figure A.7: A conductor closed by a constriction with N 
channels (quantum point contact, QPC) and with an ab¬ 
sorbing probe. 


The idea is then to estimate G by noticing that a 
conductor with one contact and an absorbing probe 


is equivalent to a two terminal conductor (Fig. A.8). 


In this latter case, the conductance is dominated 
by the resistance of the constrictions, hence I/G = 
1/(A^T) + l/{NaTa), where Na and Ta characterise 
the second contact. We deduce 


N-G = N 


NT+NaTa{l-T) 

NT+NaTa 


(A.3) 


Expansion (A.2l in powers of E and expansion 


of (A.3l in powers of NaTa can be identified by 


setting NaTa = c (E/A), where c is a dimensionless 


The relation between the dynamical problem (real a;) 
and the static problem with absorption (complex uj) was used 
in another context in Refs. |106l 11541 1141 FTTl I80| . 



Figure A.8: A two terminal conductor. 


constant. Hence 

{Tt {S\s - A)S{e + A)}) = nNNNNTa. 


T 


= iV-c- + — (C-) +G(r3). (A.4) 


Comparison with (|A.2|) and the fact that 

luce 

(A.5) 


(Tr I Q}) = 27r/A shows that c = 27r. We deduce 

8^2 


(Tr{Q2}) = (^^2):^ 


NTA^ 


leading to (113). 


Appendix B. Partial versus proper times 

Let us make few remarks on the relation between 
partial time delays and proper time delays. We con¬ 
sider a scattering situation with N channels. The 
5-matrix may be represented under the form 

S=Ue^^W (B.l) 

where 0 = diag(27yi, • • • , 2'qi^) gathers the N phase 
shifts. We deduce the following representation for 
the Wigner-Smith matrix : 

Q = G (5^0) -b i [5^ )S - UdJA^\ (B.2) 

The first term involves the partial time delays : 
5e0 = diag(fi,--- ,T 7 v), which are intrinsic char¬ 
acteristics of the scattering matrix (independent on 
the basis). The second term controlled by G, i.e. 
by the choice of basis, is the origin of the difference 
between partial and proper times. 

We now say a little bit more on the unitary ma¬ 
trix U in the N = 2 channel case. Going back to 
the notations introduced in § |2.2[ we can rewrite 
the relation between partial waves ([^ (eigenstates 
of the 5-matrix) and the left/right scattering states 
in the matricial form 

Ae,l(a;)\ _ ARAf^s,L{^)\ moN 

- (Ai.2 An,2) 
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which involves the transpose of the matrix which 
diagonalises the 5-matrix 



therefore, when eigenstates are properly normalised 
(associated to a measure de), we can write 
{(t>e,a\'fpe\b) = %a(£)d(e - s') where a € { 1 , 2 } 
and b G {L, R}. 


Appendix B.l. Symmetric case 

We first consider the symmetric case r = r' 
and t = t' [in ID this corresponds to a symmet¬ 
ric potential V{—x) = V{x)\. The left and right 
scattering states are simply related by ipe,Ri,x) = 
4 ’e,L{—x) and the two partial waves are the 
symmetric/antisymmetric eigenstates 4 >e,ai—x) = 
hence the unitary matrix is inde¬ 
pendent on the energy in this case 


U = 


1 /I 1 


72 VI -1 


(B.5) 


and the partial and proper times coincide 


Ti = fi and T 2 = f 2 . (B. 6 ) 



Figure B.9: Proper time delays (continous blue line) and 
partial time delays (dashed red line) as a function of k = y/e 
in a particular example of ID Hamiltonian with potential 
V{x) = Ai 5 (x) + A 2 5{x — a) with Ai = 2 and X 2 = 4^ (a = 1). 


Appendix B.2. Asymmetric case 

In general the matrix U carries some energy de¬ 
pendence. In order to illustrate the difference be¬ 
tween the two sets of characteristic times, we con¬ 
sider a specific ID example : we compute these 
times for the Hamiltonian H = —d^ + Ai S{x) + 
X 2 S{x — a). We consider the basis of left and right 


scattering states in which the 5-matrix takes the 
form ([^. The three coefficients r, t and r' can 
be computed easily thanks to a transfer matrix ap¬ 
proach. Diagonalisation of the 5-matrix provides 
the phase shifts and hence the partial time delays fi 
and f 2 . We also compute the Wigner-Smith matrix 
in the basis of left/right scattering states (cf. § 2.21 
and diagonalise it, we get the proper time delays ri 
and T 2 - We check that for the symmetric potential 
(Ai = A 2 ) the two sets exactly coincide, Eq. (B. 6 ). 
In the asymmetric case (Ai 7 ^ A 2 ) the four times 
are plotted in Fig. B.9 as a function of ^/e. We see 
that proper and partial times are very close (they 
get closer as the energy grows). At the crossing 
points where the two partial times are equal we can 
however observe that partial times and proper times 
differ. This follows from the fact that the second 
term of Eq. (B.2) can be understood as a perturba¬ 
tion which produces anticrossing of the eigenvalues 
of Q. 
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